Growth of three-dimensional structures by atomic deposition on surfaces containing 

defects : simulations and theory 



Pablo Jensen (a), Hernan Larralde (b), Muriel Meunier (a) and Alberto Pimpinelli (c) 
(a) Departement de Physique des Materiaux, Universite Claude Bernard Lyon-1, 69622 Villeurbanne Cedex, FRANCE; 
(b) Instituto de Fisica, Lab. de Cuernavaca, Apdo. Postal 48-3, CP. 62251, Cuernavaca, Mor., MEXICO; 
(c) LASMEA, Universite Blaise Pascal, Les Cezeaux,F-63177 Aubiere Cedex, FRANCE 

We perform a comprehensive study of the formation of three dimensional (pyramidal) structures 
in a large range of conditions, including the possible evaporation of adatoms from the surface and 
the presence of surface defects. We compare our computer simulations to theoretical calculations of 
the growth and find good agreement between them. This work clarifies precedent studies of three 
dimensional growth and predicts the island size distributions obtained in the different regimes. 
Finally, we show how our analysis can be used to interpret experimental data. 
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How can one grow useful thin films or nanostructures from atomic beams? The usual and most effective way is 
certainly by a carefully controlled and inspired trial and error method. Theoreticians dream that another possibility 
may exist : by knowing the detailed atomic mechanisms that govern thin film growth, one could in principle tailor the 
morphologies to the desired application. Fortunately, there are many other justifications to the study of these atomic 
mechanisms : for example the understanding in terms of atomic mechanisms of growth experiments carried under 
controlled conditions, where great care is taken to avoid complications (contamination, surface defects, etc. see for 
example |l - 6) ) . To be able to interpret more complex situations demands the study of models including many atomic 
processes [ 7 1 . The first models only included deposition, diffusion of the adatoms and their irreversible aggregation 
to form flat islands |8|-|l6|. Ratsch et al. |l7j] improved these models by including reversible aggregation in order to 
reproduce the formation of compact islands. In this paper we get closer to new experimental situations by including 
two new ingredients : the growth of three-dimensional islands and the presence of surface defects which act as traps 
for the monomers, in presence of adatom evaporation. 

Evaporation, i.e. the possibility of desorption of adatoms from the surface, is a feature that should be observed 
for any system at high enough temperatures. In this sense, it is a phenomenon that is as general as the rest of the 
ingredients of recent models of film growth, and, as we have shown in the case of two dimensional (2D) growth [fj"8| , 
is capable of completely changing t he qua ntitative behaviour of the system. Moreover, evaporation is present in a 
number of experimental situations jl 19 2(J, where 3D islands are built. We should add that thin film growth models 
which include evaporation have already been studied using a mathematical analysis of rate equations ]Tc| , p^ - |T6| , pO| -p3| . 
Computer simulations of such models have, to our knowledge, never been carried out. The point is that computer 
simulations represent an "exact" way of reproducing the growth, in the sense that they avoid the mean-field approx- 
imations of rate-equations approaches Jl3|,p4|,|5| . We have shown previously that the mean- field equations could 
lead to wrong predictions in the case of 2D growth : is this also the case in 3D growth? 

Three-dimensional (3D) structures are often observed in the growth of heteroepitaxial and non-epitaxial films. A 
simple explanation of the formation of such structures, based in the bonding energies between adsorbed atoms and 
adsorbed atoms and the substrate can be found in |^6| . We note that such thermodynamic arguments are not always 
trustworthy, since kinetics play an essential role in determining the growth morphologies |2j]-|29j . We will not consider 
this point here and will just simulate 3D growth in a schematic way, forcing the adatoms to build pyramidal islands 
as they aggregate on the surface. A more realistic model should reproduce the geometric structure of the islands as a 
result of the different relevant energies (edge diffusion around the island, barrier for down and up-hill diffusion, . . . ). 
This is beyond the scope of the present paper where we only want to investigate the consequences of the formation of 
3D structures on the growth characteristics (mainly saturation island density and island size distribution as a function 
of the growth parameters) . The possible influence of surface defects has been stressed repeatedly [^,|o[^l) . We want 
to clarify the manifestations of defects on the growth and to check simple mathematical analysis of the growth by 
computer simulations. 

We study here the first stages of the growth, roughly until the number of islands on the substrate saturates. The 
reason is that it is in this stage that such a simplified model can be of some help to experimentalists who want to 
understand the microscopic processes present in their experiments. These fundamental microscopic processes are most 
easily detected in the first stages of the growth, since in the subsequent stages additional processes can be involved 
(additional diffusion paths, interlayer transport, geometrical details of the lattice . . . ). 

The paper is organized as follows. Section I briefly presents the model and discusses some of its approximations. 
Then, in section II, we study the growth of 3D islands, first by a simple scaling analysis in the absence of defects, then 
by a more rigorous mathematical mean-field treatment, where we also include the influence of surface defects which 
act as perfect traps on the surface. In section III, computer simulations are used to test these analytical predictions 
and to calculate the island size distributions in the different cases. After a discussion (section IV) where we compare 
our analysis to precedent studies, we show in section V how experimental results can be analyzed using these results. 



I. PRESENTATION OF THE MODEL 



In this work we will describe the properties of a still oversimplified submonolayer thin film growth model which 
includes five important physical ingredients for these systems: 

(1) Deposition. We will assume that atoms are deposited at randomly-chosen positions of the surface at a flux F per 
unit surface per unit time. Atoms that happen to fall on the islands that are formed on the surface are accomodated 
at their proper pyramidal position (see below). 

(2) Diffusion. Isolated adatoms can move in a random direction by one diameter, or one lattice spacing, which 
we will take as our unit length. We denote by r the characteristic time between diffusion steps and D = 1/ (4t) the 
diffusion coefficient (the atomic size is taken as the length unit). 
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(3) Evaporation. Isolated adatoms can evaporate off the surface at a constant rate. We denote by r e the mean 
lifetime of a free adatom on the surface. It is also useful to define X$ — y/Dr e , the mean adatom diffusion length 
before desorption 

(4) Aggregation. If two adatoms come to occupy neighboring sites, they stick irreversibly and form an island. As 
more adatoms are captured, the island develops as a pyramid (see below). Islands are assumed to be immobile and 
do not evaporate. 

(5) Defect trapping In some parts of the paper, we introduce a concentration c of "defects" on the surface. These 
defects, which are randomly distributed on the surface, act as perfect traps for the monomers. Therefore, a monomer 
which happens to occupy a defect remains there forever and serves as a nucleation center for island growth. 

In the following, we call particles or adatoms the isolated atoms (or monomers) that are deposited on the surface, 
and islands a set of connected particles (thus excluding the monomers). 

Some remarks on the assumptions of this simple model regarding its connection to the experiments are now ad- 
dressed. 

(a) Island morphology — We force the islands to assume a pyramidal shape by immediately moving an adatom that 
touches an island to the desired position (Fig. ^). The pyramidal shape is adopted because we wish to understand 
experiments carried out by one of us (MM) on the system Pd/NaCl Q where the islands take approximately such a 
shape. We note that this assumption does not affect crucially the growth : it should not affect the growth exponents 
which are only determined by the fact that islands are three-dimensional (i.e. their mass increases as their radius to 
the third power, see below). 

(b) Island diffusion — We neglect in this model the possibility for dimers, trimers or larger islands to diffuse on the 
substrate. Island diffusion has been observed in experiments p2] and molecular-dynamics simulations [33|. The effects 
of 2D island diffusion on the growth of thin films without evaporation have been addressed in Refs. ]i"2|, 24 ,p5p^p^,|35[| . 



II. MATHEMATICAL ANALYSIS OF THE GROWTH OF 3D ISLANDS 



A. Qualitative description 

Before going into the details of the calculations and their confirmation by computer simulations, we present a simple 
picture of the growth mechanisms of the submonolayer structures under consideration. We are interested mainly on 
two points : the time evolution of the island concentration on the substrate and the island concentration at saturation 
as a function of the growth parameters : flux, diffusion and evaporation times and defect concentration. 

The qualitative evolution of the system is essentially as follows. The system initiates as a clean empty surface. 
Monomers are then deposited at a constant rate on the surface and are allowed to diffuse and/or desorb (evaporate). 
When two monomers meet, they aggregate irreversibly to form a static island (an island is also created when a 
monomer is trapped by a defect) : this is island nucleation. As more of these encounters occur, the island density 
increases with time. Competing with this nucleation process, monomers are also captured by islands which then 
become larger. At some point, islands are so large that they quickly capture the free monomers, which reduces the 
chances of two monomers meeting to nucleate a new island. Therefore, the number of islands grows less rapidly. 
Moreover, when islands become large, they start touching (coalescing). These two effects lead to a saturation in 
the number of islands. Interestingly, the saturation is attained when the surface coverage reaches a value close to 
.15, independently of the parameter values. This is equivalent to saying that saturation is reached when the mean 
island radius R is a fixed fraction of the island- island distance (the coverage 9 is given by 8 ~ NR 2 ~ (R/la) 2 where 
la is the mean island-island distance). A more detailed discussion of the evolution of the systems can be done by 
distinguishing several cases according to the relative importance of the different phenomena : diffusion, evaporation, 
defect concentration. It is useful to define a typical length scale for each of these processes: lec = {Ft) 1 / 7 is a typical 
distance between islands when evaporation and defect concentration are not important (see below), Xs = \JD>T e is, 
as defined above, the adatom diffusion length before desorption and t& e f ~ 1/c 1 / 2 is the typical distance between 
defects. Now, depending on the relative values of these three length scales, growth will be dominated by different 
mechanisms which we turn on to describe in more detail. Note that lec corresponds to the typical island-island 
distance at the saturation time, i.e. it is not the actual island-island distance (for example at the beginning of film 
growth). Therefore our following qualitative discussion is only approximated, and the more precise calculations of 
later sections are necessary to justify it. 
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1. "Dirty" substrates : high defect concentration 



By "dirty" substrate, we mean that the island concentration is affected by the defect concentration. We will show 
that this is true only if ice is much larger than idef, he. the defect concentration is high enough, even in the absence 
of evaporation (if evaporation is present, it can only decrease the number of islands, therefore increasing the relative 
importance of defect concentration). The simplest case is when the defect concentration is very high, namely if idef 
is much smaller than both X$ and tec- I n this case, island nucleation is completely dominated by the trapping 
of adatoms by defects, which leads to two main effects : first, the island concentration reaches its saturation value 
rapidly (roughly in a time c/F); second, this saturation value is AT sat = c, i.e. all defects are occupied by islands but 
there are no islands created elsewhere. 

Note that this situation corresponds to a "low" temperature case, when X$ is large enough. At higher temperatures, 
one could have Xs -C idef (but still idef <C tec to remain in the "dirty" substrate case). In this case, the monomer 
concentration is dominated by evaporation and island nucleation still occurs on defects. The saturation density is 
still equal to the defect concentration but the kinetics is different : the time needed to reach saturation is roughly 
f c ~l/(F(l + Xf)). 



2. "Clean" substrates : low defect concentration 



Here we study the cases for which the island concentration is not affected by the presence of defects, when the 
substrate is "clean" enough. This is true when ice is much smaller than idef, irrespective of the Xs value. In principle, 
three cases can then arise, depending on the relative value of Xs against these two lengths. Our calculations will 
show (section [IB) that there are only two asymptotic regimes : complete condensation and high evaporation. In the 
complete condensation (CC) case, adatoms do not evaporate from the surface and island growth proceeds mainly by 
capture of the diffusing adatoms. On the contrary, in the high evaporation limit, the number of adatoms is limited 
by evaporation and the most important growth mechanism for islands is that of direct impingement of adatoms on 
top of them (the contribution from the adatoms diffusing on the surface is negligible). 

The first case corresponds to ice <C Xs, independently of the relative ordering of Idef and Xs- Then, adatoms 
almost never evaporate before aggregating (after this, they are safe since islands do not evaporate). The monomer 
density rapidly grows, leading to a rapid increase of island density by monomer- monomer encounter on the surface. 
This goes on until the islands occupy a significant fraction of the surface, roughly 1%. Then, islands capture rapidly 
the monomers, whose density decreases. As a consequence, it becomes less probable to create more islands, and we 
see that their number increases more slowly. When the coverage reaches a value close to 15%, coalescence will start 
to decrease the number of islands. The maximum number of islands N sat is thus reached for coverages around 15%. 

The second case corresponds to the opposite situation : ice 3> Xs, with still ice ^ &def- This happens when r e 
is small, for example at high temperatures. In this regime, evaporation significantly alters the growth dynamics. The 
main point is that now the monomer density becomes roughly a constant (Fr e ), since it is now mainly determined 
by the balancing of deposition and evaporation. Then, the number of islands increases linearly with time (the island 
creation rate is roughly proportional to the square monomer concentration) . We also notice that only a small fraction 
(1/100) of the monomers do effectively remain on the substrate, as shown by the low sticking coefficient value at early 
times (the sticking coefficient is the ratio of particles on the substrate (the coverage) over the the total number of 
particles sent on the surface (Ft)). This can be understood by noting that an island of radius R grows by capturing 
only the monomers that are deposited within its " capture zone" (the circle of radius R + Xs centered on island's 
center). The other monomers evaporate before reaching the islands. As in the case of complete condensation, when 
the islands occupy a significant fraction of the surface, they capture rapidly the monomers. This has two effects : the 
monomer density starts to decrease, and the sticking coefficient starts to increase. Shortly after, the island density 
saturates and starts to decrease because of island-island coalescence. 

Note that one could have in principle Xs < 1, i.e. the adatoms evaporate before diffusing. This situation, although 
apparently uncommon, is not physically impossible and it also allows us to test our predictions over a larger range of 
parameters. We note that, in contrast to what is observed for two dimensional (2d) islands [Q, particles deposited 
on top of islands significantly contribute to island growth. In presence of strong evaporation, this is actually the 
dominant mechanism for island growth, whereas for the 2d case this happens only in some special regimes 18 . 
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3. Summary of our results 



We present her e th e su mma ry of our results in the different limits described above. These results are derived in 



detail in sections [I B and [I C . For each regime, we give in the order the saturation island density N sat , the thickness 
at saturation e sa t (i.e. the thickness when the island density first reaches its saturation value), the thickness at 
coalescence e c (i.e. the thickness when the island density starts to decrease due to island- island coalescence), and the 
scaling kinetics of the mean radius as a function of time before the saturation island density is reached. We recall 
that lec = (Ft) 1 / 7 and X s = y/rjr'. 

Clean substrate 

high evaporation : Xs <C lec "C Idef 
N sat ~ [Fr e (l + X 2 )] 2 / 3 
e sa t~e c ~ [Fr e {l + X 2 s )}- 1 ^ 
i?~ Ft 

low evaporation : lec -C Xs <C Idef or l C c <C Idef -C Xs 
N ~ (£) V7 

1 y max \ D J 

/mi/7 

&sat ^ C c ^ \ p ) 

R ~ (FDt 2 ) 1 / 9 ~ t 2 / 9 



Dirty substrate 

high evaporation : Xs <C Idef "C lec 
e sat ~ rnrwi 



e c 

R^ Ft 



l 

1/2 



low evaporation : i^e/ "C Xs <C Zee or Z de / <C Zee -C -X's 

Nmax ^ C 
dsat ~ C 

e l 

C C c l/2 

i? ~ (— ) for t < c/F, i.e. before saturation 



R ~ (Ft/c) 1 ' 3 between saturation and coalescence (c/F < t < 1/Fc 1 / 2 ), see section [I C 



B. Scaling arguments for defect- free ("clean") surfaces 

In this section we present simple scaling arguments that allow to find the dependence of the maximum island density 
N m ax as a function of the deposition parameters (Flux F, Diffusion time r and Evaporation time r e ). These arguments 
were originally formulated in jl0| for the special case of growth of 2D islands by atom deposition on a high-symmetry 
terrace, neglecting evaporation. Here, the argument is extended to the case of non-negligible evaporation. We recall 
that the atomic size is taken as the length unit. 

The first stage of the argument requires the determination of the nucleation rate per unit surface and time, l/r nuc . 
A nucleation event takes place when an adatom meets a critical island of size s*. We call N s * the density of critical 
nuclei (clusters of size s*) and, following Refs. [ p2p5) ], we assume that N s * satisfies Walton's relation N s * ~ p s . 
Thus, 

— ^(F + Dp)N s , (2.1) 

7~nuc 
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where D — 1/(4t) is the adatom diffusion constant. The term FN S * accounts for direct impingement of atoms onto 
critical islands and the second term for critical island growth by monomer diffusive attachcmcnt. 

Another, independent equation can be written down to relate the nucleation rate and the stable (s < s*) island 
density N. It states that in the area l 2 s — 1/N occupied by an island, only one (on average) nucleation event takes 
place, during the time t c needed for the growing islands to come into contact. Thus, 

(2.2) 

7"nuc 

The time t c is readily computed by knowing the growth velocity of an island, which in turn requires the knowledge 
of the adatom density. We consider in the following three situations of interest for this paper. 



1. Negligible evaporation 



The adatom density results in this case from a balance between deposition at a rate F and capture by the stable 
islands at a rate DpN, so that 



p ps F/(DN) « F£ 2 JD 



(2.3) 



The rate of growth of the volume of an island of linear size R is diffusion-limited in this case, so that d(R 3 )/dt ps Dp 
and R 3 ps Ft/N. At t = t c , R ps £ s , and thus 



t c ps Nl\/Fm 1/{FN 1/2 ) 



(2.4) 



Using (2.1) and (2.3), one finds 



D[F/(DN)} 



s* + l 



(2.5) 



(here, th e di rect i mpin gement term is negligible). 
From dJh and (O), 



— « FN 3 / 2 

''nuc 



Finally, (2.5) and (2.6) yield [CT2J 



JV : 



2s*/(2s*+5) 



(2.6) 



(2.7) 



2. Strong evaporation 



Strong evaporation means the adatoms are more likely to disappear due to desorption than to be captured by 
an island. In other words, the adatom diffusion length before desorption, Xs = v 'Dr e , is shorter than the average 
island- island distance, t s . In this case, the adatom density results from a balance between deposition and desorption 
at a rate p/r e , so that 

p^Fr e ^ FX 2 S /D . (2.8) 

A 3-D island grows by two mechanisms in the case of strong evaporation: (i) by capture of the adatoms falling on 
the surface at a distance smaller than Xs; (ii) by direct capture of all adatoms falling on it. Thus, dR 3 /At ps F(XsR+ 
X 2 + R 2 ) ps F{R 2 + X S R). Solving for R with the condition R = at t = 0, one gets R-X s ln(l + R/X s ) = Ft. 
At t = t c , R ~ i s > Xs, which means that direct capture always dominates. It follows Ft c ps £ s , or 

t c Kl/[FN 1/2 ] . (2.9) 
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Using (2.1) and (2.8), one finds 



From (2.2) and (2.9) 



« (F + Dp)p s ' a F(l + X 2 )p s * a F(l + X 2 )(FT e ) s 

'mxc 



(2.10) 



(2.11) 



Finally, ( |2.10[ ) and ( p.lip yield 

A^a (Fr e ) 2s * /3 (1 + X|) 2/3 
We discuss the preceding results in section [j^. 



(2.12) 



3. Crossover between the two preceding regimes 



We will show in this section that it is possible to derive interpolation formulae describing the crossovers between 
the no-evaporation and the strong-evaporation regimes. 

To do this, we will use Burton, Cabrera and Frank's theory of step flow j36|. The adatom density will be computed 
on a terrace bounded by two parallel steps, at a distance i. They may be the steps in the orderd array of a vicinal 
surface; in this case, we will let I = d. Or they may represent the edges of a big enough island; in the latter case we 
will let i = t s = 1/VN. 

The adatom density obeys the equation 



F + DW 2 n 



(2.13) 



The time r is the average lifetime of an adatom. Since adatoms disappear either by capture by islands, or by 
desorption, we will alternatively let 1/t = DN = D/l 2 s or 1/t = l/r e = D/Xg. In both cases, the notation 
k 2 = 1/(Dt) will be used. 

In the quasi-stationary approximation pq|, h a 0, and choosing the origin x — at the terrace centre, equation 
( [2.13 ) can be solved with the conditions n{±£/2) — at boundary steps. The solution reads 



n(x) = Ft 



1 



cosh(Kcc) 
cosh(«^/2) 



(2.14) 



This formula will be nee ded t o compute the nucleation rate (2.1). The latter is an average quantity, independent 
of x. We let thus x = in (2.14), since the higher nucleation probability is at the terrace centre, given the symmetry 
of our problem. One finds 



n = Ft 



1 - 



1 



cosh(«^/2) 



= 2Ft 



sinh 2 (/^/4) 
cos1i(k£/2) 



(2.15) 



where we used the identity cosh(x) — 1 = 2 sinh 2 (a;/2). 

The next task is the determination of the island density. To this end, it suffices to consider the total island density 
N. Its time variation is simple: N increases each time a new island is nucleated, so (N)i = DpN s * . On the other hand, 
N decreases when two islands touch and coalesce. Following the authorities J2^,|2j|, we write (N)2 = —(dA/dt)N 2 , 
where A ~ R 2 is the average area of an island of linear size R. This means that coalescence results from binary 
encounters of immobile islands, whose area increases at a rate dA/dt. Collecting (N)i and (N)2 yields 



N = DpN s « - ^N 2 
dt 



(2.16) 



At stationarity, which is what we are interested in, N = 0. Note that, by definition, coalescence takes place when 
A ~ £ 2 . One can thus write 



DpN s 



dA\ 
dt ) 



N 



(2.17) 
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The final task concerns therefore the evaluation of the growth rate of an island. This can be done by noting that 
the mass M. of an island increases either by direct capture of atoms from the beam, or by diffusion of adatoms on the 
surface. Since the surface diffusion current of adatoms is — DVn, assuming circular (d = 2) or hemispherical (d = 3) 
islands of radius R one can write 



AM 
~df 



FR 



d-l 



D 



dn 
dr 



(2.18) 



The result can be found in a number of papers [|18|,|22 23 39]. It makes use of Bessel functions, which are not easy 
to manipulate. Approximate results, which have the merit of reproducing the correct limiting regimes (see below), 
and of allowing analytical expressions to be writte n down, will be used. 

The adatom density is still given by Eq.( 2.15| ), and the nucleation rate by Eq.(2.1). The growth rate of the 
(projected) area of a 3-D island has two contributions, as for 2-D islands: a diffusion-limited contribution, which is 
given by 



dA 
dt 1 



FXs4tanh(«4/2) ; 



and a direct capture contribution, 



dA, 
~dt 



\R=t, 



FL 



A useful interpolation formula between ( 2.19| ) and ( 2.2C ) is 



dA, 



-^Ir^s « Ft. ( 1 + it ) tanh( K 4/2) . 



Finally, (|]| and ( fjj| ) yield 

(i + XsY 



Ft, 



smh 2 (ne s /4) 



cosh(re£ s /2) 



s*+l 



F<„ [ 1 + y~ ) tanh( K 4/2)A 2 



or, 



(l + X s ) 2 (Fr e Y 



cosh( K A- 1 / 2 /2) 



sinh 2 ( K A- 1 /2/4) 



s*+l 



tanh(fdV~ 1/2 /2)(l + X S N 1/2 )N 3/2 



Multiplying both sides by A| yields 



Xf +3 (1 + X 2 ){F/D) S * » fiXgN 1 / 2 ) , 



where 



Letting A 3D cvap = {FX%/ 'D) 2s ' 73 (l + Xf)" 1 " , formula (|2_2j) can be cast in the form 

f3Dcv&p(XgN) = XgN3j) eV3l p , 

where /3Devap(^) = [fix 1 / 2 )} 2 / 3 . Inverting f(x) yields the crossover scaling function 

XgN = 53Dovap(A^iV3Dcvap) • 

The function g3Dcvap(aO has the following properties: 

f x 3/(5+2s*) f Qr 3. 



cosh(^) 



~1 s* + l 



sinh 



2x , 
2 I 1 



tanh ( 2~ ) (1 + x)x 3 



2/3 



53Devap(x) 



for x — > 0. 



(2.19) 
(2.20) 
(2.21) 

(2.22) 

(2.23) 
(2.24) 
(2.25) 

(2.26) 
(2.27) 

(2.28) 
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C. Rate equations 



In this section we study the growth in presence of defects using rate equations [B7| . 

We will consider that the islands are semi-spherical "droplets" as a result of "fast" large scale reacomodation of 
the monomers forming the island. In addition to evaporation, we will also consider the effects due to the presence 
of point defects on the surface. We will assume that these defects serve as perfect nucleation points, and so, as the 
concentration of defects increases, the system passes from homogeneous to heterogeneous nucleation. 

We can write the evolution of the density p of monomers on the surface as: 

^=F-p/r e -Fp- a a {2p + c def ) - a t N. (1) 

The first term on the RHS denotes the flux of monomers onto the surface. The second term represents the effect 
of evaporation. The third term is due to the possibility of losing monomers by effect of direct impingement of a 
deposited monomers right beside a monomer still on the surface to form an island. The next terms represent the loss 
of monomers by aggregation with other monomers, nucleation on defects and aggregation with islands respectively. 
The factors a and Ci are the "cross sections" for encounters and are detailed below. 

The number iV of islands will be given by: 

dN 

— = F{p + c d ef) + cr (p + c def ) (2) 

where the first term represents the formation of islands due to direct impingement of deposited monomers onto other 
monomers or defects, and the second term accounts for the formation of islands by the encounter of two monomers 
diffusing on the surface, or a monomer on a defect. It should be noted that since this description yields scaling results, 
some numerical factors are omitted. 



The concentration of "free defects" varies as 

dCdef 

dt 

The total island mass density (M) changes as 



-(a + F)c de f- 



^L = 2[Fp + a oP ] + a t N + a c def + FR 2 N, (5) 



where the direct impingement term is proportional to the area of the islands. 
The typical island radius R will be given by 



/M\ 1/3 

r ~[n) ' ^ 

and the actual island coverage will then be 

Q ~ NR 2 = (NM 2 ) 1/3 . (7) 
The expressions for the diffusive cross sections in the presence of evaporation were calculated in |H| . 

2g for R»X S , 

(12) 

Dp forR«X s . 

Where X$ = y/Dr e is the typical distance traveled by a diffusing particle on the surface before it desorbs. And: 



<Jo ~ Dp. (13) 

These mean-field cross sections only depend on the radius of the sphere touching the surface, and not on the height 
of the island, therefore they are the same for 2d or 3d particles of same radius. 







As the most we can expect to obtain from this description are the scaling behaviors, we focus only on the extreme 
limiting cases of the system described by the rate equations. 

For the case of negligible evaporation and high enough initial concentration of defects, we expect the equations to 
reduce to 

±~ F -Dp(c def + N). 
dN 

— <~ Dp(p + c def ) 



Adding the second and last equations, we have: 



dCdef „ 

— = -Dpc def . 



d(N + c def )_ Dp2 



dt 

For high enough initial concentration of defects c, we can write 



N(t) ~ c - c def (t) + D(-j^) 2 t ~ c 



P ~Wc 



Cdef(t) ~ ce 

The island density at saturation is therefore equal to the initial defect concentration c. To find what is "high 
enough" for the defect concentration, we note that nucleation events start to reduce the number of islands roughly 
when the island size is of the order of the distance between islands, or equivalcntly, when the coverage reaches a 
constant value. Since the deposited mass grows as Ft (evaporation is negligible), the surface coverage will be given 

by 

fi ~ N^iFtf/* ~ (cF 2 t 2 ) 1/3 
Thus, the coalescence time is t c <~ l/(Fc 1 / 2 ) and the cross over to the "clean system" occurs when 

(^)~c i.e. c~{F/D)W. 

It should be stressed that in this defect dominated, or "dirty", regime the island number density saturates far before 
coalescence, indeed the saturation time t s can be estimated by the characteristic time for defect ocupation: t s ~ c/F. 
After t s , the island density remains essentially constant and the islands grow (the typical island radius can be easily 
shown to grow as R <~ (Ft/c) 1 ' 3 ) until coalescence. 

At high evaporation rates the aggregation of mass on the surface is dominated by "direct impingement" of particles 
on the islands, and the concentration of monomers on the surface is determined by the balance between deposition 
and evaporation. Under these conditions, at high initial defect concentration we will have: 

p~Fr e 

^ - FT e Dc def + F 2 r e + DF 2 r 2 e 
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dc def (t) 

— -j- F(Dr e + l)c def 

These equations can be solved immediately, from which we get: 

and 

N ~ c[l - e-I^+H^ *]] + F 2 r e [X 2 s + l]t, 

where we have substituted X s = \fDr e . 

From these expressions it is evident that the time at which N saturates is t s ~ 1/(F[1 + X 2 }), after which time, the 
island density reaches the value N ~ c (this, of course, assuming that the homogeneous contributions are negligible 
throughout the evolution of the system). It is then straight forward to find that the mass deposited on the system 
grows as M ~ cF 3 t 3 and the typical island radius R ~ Ft. The coverage increases as ft ~ cF 2 t 2 , from which we can 
estimate the coalescence time t c ~ l/Fc 1 ^ 2 . While it can be imagined that these times could occur in the wrong order, 
we have checked that for this to be the case, an initial concentration of defects larger than one would be required. 

By comparing the maximum value of the subdominant homogeneous term with c, we can determine that this high 
evaporation dirty regime is attained when c >> (Ft e [1 + X 2 ]) 2 / 3, . 

Finally, for completeness, we sketch the derivation of the results for the "clean" substrate with negligible evaporation. 
Under these circumstances, the monomer density on the surface is determined by the balance between deposition and 
diffusive capture by the islands on the surface. Thus, the monomer and island densities are given by 

F , dN 2 F 2 

and — ~ Dp" - 



' DN dt r DN 2 ' 

From these equations we find that the island density is N ~ (F 2 t/D) 1 ' 3 . Also, as we are considering the case of 
negligible evaporation, the mass deposited on the surface will be M ~ Ft. From these quantities we can infer the 
behavior of the typical island radius to be given by R ~ (M/iV) 1 / 3 ~ (DFt 2 ) 1 / 9 , and the saturation and coalescence 
times to be t c ~ t s ~ (D/F 8 ) 1 / 7 , at which times the maximum value of N is N max ~ (F jD) 2 ! 7 , as obtained in the 
previous section. 



III. COMPUTER SIMULATIONS 



In the following paragraphs, we test the assumptions and predictions of the analysis given in the preceding sections. 
We also show results that are not attainable from this mean-field calculations, namely the island size distributions. 

Our computer simulations generate sub-monolayer structures using the four processes included in our model (see 
the introduction). Here we take r = 1 as the time scale of our problem. The monomer diffusion coefficient is then 
given by D — 1/4. We use triangular lattices (six directions for diffusion) of sizes up to 2000 x 2000 with periodic 
boundary conditions to limit finite size effects. 

The program actually consists of a repeated loop. At each loop, we calculate two quantities Pdrop = F/(F+p(^ + ^)) 
and pdif = (p/t)/(F + p(— + -)) that give the respective probabilities of the three different processes which could 
happen : depositing a particle (deposition), moving a particle (diffusion) or removing a particle from the surface 
(evaporation). More precisely, at each loop we throw a random number p (0 < p < 1) and compare it to pdrop and 
Pdif - If p < Pdrop, we deposit a particle; if p > p dr0 p+Pdif, we remove a monomer, otherwise we just move a randomly 
chosen monomer. After each of these possibilities, we check whether an aggregation has taken place and go to the 
next loop (for more details, see [0)- 



1. Checking the crossover scaling 



Before looking in detail into the different regimes predicted in section II C, we summarize our simulation results in 
Fig. |[ We show there all our data for N max as a func tion of the parameters. Our scaling analysis predicts that the 
data should fall into a single curve, given by Equation 2.28| . We see that the data remarkably confirms our analysis, 
over more than 30 orders of magnitude. This gives us confidence on our entire approach and its predicted exponents, 
which we now turn on to check in more detail. 



We now check that the results summarized in section [I A 3 arc correct 
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2. Scaling of the maximum island density as a function of incident flux 



Figure y shows the evolution of the maximum island density as a function of the flux for different evaporation 
times. Each of these curves is different from the others, since they correspond to different evaporation times. However, 
according to our preceding analysis, they should all present a transition from the low evaporation regime to the high 
evaporation regime. This can be detected by a change of slope, from N max ~ F 2 ! 3 in the high evaporation regime 
(solid line) to N max ~ F 2 / 7 in the low one (dashed line). Of course, this regime change does not occur for all the 
curves at the same value of the Flux, since the para meter that determines that change is not the Flux but rather 



Xg = t c /t. Figure |3| shows that the results of section II A £ accurately describes the behaviour of our model, at least 



concerning the Flux evolution of the maximum island density. We now turn to the other variable, the evaporation 
time. 



3. Maximum island density as a function of evaporation time 



We show in Figure |j the dependence of the maximum island density on r e (t = 1). We notice that for high 
enough evaporation times, the island density tends to become r oughly constant, as predicted by our calculations. For 
lower values of r e , N max changes rapidly. We predict (section QA3) two regimes : for 1 <C t £ -C F -1 / 3 , we expect 

Nmax ~ Te^ 3 , while for T e <C 1, we expect N max ~ r e 2 / 3 . This last regime is clearly seen for the curves obtained for 
fluxes F = 10~ 6 and F — 1CP 4 (squares and diamonds respectively, the slope 2/3 is given by the solid line). The first 
regime is difficult to see for two reasons. First, we need X$ ^> 1 and strong evaporation, i.e. I ^> Xs- This means a 
very low island density, meaning very long computing times and large lattices. Second, the crossovers with the two 
other regimes (exponents 2/3 and 0) tend to mask the exponent 4/3. Taking a lower value for the flux (F — 10 -8 , 
filled circles), we can see that the slope in this intermediate regime is larger than 2/3. 



4- Mean island radius versus time for clean substrates 

Our treatment predicts two limiting regimes for the power-law (r ~ t°) evolution of the mean island radius as a 
function of time : /3 = 1 in cases of strong evaporation and [3 = 2/9 for complete condensation (we only treat here 
the case of "clean" substrates"). Fig. ||a shows that we observe indeed an exponent very close to 2/9=0.22 when 
evaporation is negligible (Xs = 10 9 , Ice — 37, Xs ^> Ice), while the exponent is close to 1 in the opposite case 
(Xs = 1, Ice = 14, X s <C Ice), see Fig. |b. Of course, intermediate cases can arise in experiments and the exponent 
is between these two extreme values, with values around 0.5-0.6 as shown in Fig. ||c (Xs = 10, lec — 27). Note that 
we have defined here the radius as (M/N) 1 / 3 where M is the total mass present on the substrate and N the island 
density, but we have checked that similar exponents are measured if one defines the mean radius as fl/N 1 ^ 2 where Q 
is the surface coverage. 



5. Dynamical evolution of island density 



Here we investigate how the different microscopic mechanisms can affect the growth kinetics. This can be an 
important help for experimentalists seeking information on which processes are actually present in their experiments 
p8| . Fig. | confirms our analytical analysis and shows that evaporation or the presence of surface defects can 
significantly alter the time evolution of island density. If defects are present, monomers will be trapped by them 
at the very beginning of the growth and the number of islands equates rapidly the number of defects, whatever the 
diffusivity of the atoms. If evaporation is present, the opposite effect is observed : since many atoms do not contribute 
to the growth (they evaporate before reaching an island), the saturation is reached for very high thicknesses (typically 
e sat > XML). 



6. Island size distributions 



Island size distributions have proven very useful as a tool for experimentalists to distinguish between different 
growth mechanisms (|^^^]. By size of an island, we mean the surface it occupies on the substrate. For the "three 
dimensional" particles studied here, their projected surface is the easiest quantity to measure by microscopy. Note 
that the projected surface for a given mass depends on the precise shape of the islands, which is assumed here to 
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be pyramidal (close to a half-sphere). Size distributions are normalized by the mean island size in the following way 
: one defines p(s/s m ) = n s /N t as the probability that a randomly chosen island has a surface s when the average 
surface per island is s m = 0/N tl where n s stands for the number of islands of surface s, N t is the total number of 
islands and 9 for the coverage of the surface. It has been shown that by normalizing the probabilities and plotting 
Sm *p{s/s m ) against s/s m , one obtains a "universal" size distribution independent of the coverage, the flux or the 
substrate temperature for a large range of their values. These size distributions can be obtained from the simulations 

HIIpl. 

Fig shows the evolution of the reseated island size distributions for three dimensional islands (pyramids) in 
presence of evaporation. We recall that size means here the projected surface of the island, a quantity which can 
be measured easily by electronic microscopy. We note the same trends as for 2d islands It is clear that the 

distributions are significantly affected by the evaporation, smaller islands becoming more numerous when evaporation 
increases. This trend can be qualitatively understood by noting that new islands are created continuously when 
evaporation is present, while nucleation rapidly becomes negligible in the complete condensation regime. The reason 
is that islands are created (spatially) homogeneously in the last case, because the positions of the islands are correlated 
(through monomer diffusion), leaving virtually no room for further nucleation once a small portion of the surface is 
covered (6 ~ 0.05). In the limit of strong evaporation, islands are nucleated randomly on the surface, the fluctuations 
leaving large regions of the surface uncovered. These large regions can host new islands even for relatively large 
coverages, which explains that there is a large proportion of small (s < s m ) islands in this regime. 

Fig |8| shows the evolution of the rescaled island size distributions for pyramidal islands nucleating on defects. Two 
main differences can be noted. First, the histograms are significantly narrower than in the preceding case, as had 
already been noted in experimental studies |H| . This can be understood by noting that all islands are nucleated at 
almost the same time (at the very beginning of growth). The second point is that the size distributions are sensitive 
to the actual coverage of the substrate, in contrast with previous cases. In other words, there is no perfect rescaling 
of the data obtained at different coverages, even if rescaling for different fluxes or diffusion times has been checked. 



IV. DISCUSSION 



Other authors have analyzed similar mean-field rate equations to find the growth dynamics and maximum island 
density in the presence of evaporation [20|-|23]| . They have also obtained different regimes identified roughly in the 
same way as in the present work. We think that there is some ambiguity in their definition of the different regimes. 
A regime - meaning a single relation between N max and the deposition parameters - should be defined only by the 



values of F, r e and r, as in Section [I A 3. This is what is required from the experimental side : given some values 
of the parameters, what will happen on the surface? Instead, previous works have introduced other parameters, such 
as the coalescence coverage, the capture cross sections, or even the island density itself, in the characterization of the 
regimes. 

Besides this general remark, we note that there is some disagreement about the different regimes between the 
various authors. Stoyanov and Kashchiev ]23| find two regimes which correspond to our complete condensation and 
strong evaporation (with Xs S> 1) cases. We have added here the case Xs <C 1. Venables et al. Q find three 
different regimes. The two extreme regimes coincide with Stoyanov's and ours, but their intermediate regime (which 
predicts N rnax ~ (FT e ) 2 ^ 5 is not observed in our simulations : this is particularly clear in Fig. ^ which shows no 

T,?^ 5 dependence : instead, the intermediate regime shows an exponent greater than 2/3, as predicted by Eq. 2.12). 
Actually, Venables's intermediate regime seems physically strange because there is no dependence on the diffusion 
coefficient of the adatoms. 

The time evolution of the mean radius has also been studied previously (for a recent review, see p2| and references 
therein). Here, we have shown that the exponent of the radius versus time power-law is 1 when evaporation is 
important (islands grow only by direct impingement), 1/3 when condensation is almost complete and the number 
of islands is constant (for example when their concentration reaches the defect concentration for "dirty" substrates) 
or 2/9 for clean substrates and complete condensation (here the number of islands is never really constant since 
nucleation does not stop until coalescence). 

In previous studies 42 1, it has been predicted that the radius shows a power-law dependence with an exponent 1/3 
for complete condensation growth and 1 in the case of strong evaporation. Intermediate values were found thanks to 
the numerical resolution of mean-field equations of island growth 42 1. 

The experimental values are rather in the range 0.21 - 0.30 for the Pd/MgO(100) system ]23] and around 0.3 in other 
systems. While the intermediate exponents (between 1/3 and 1) can be explained as pertaining to the intermediate 
regime (see section [II 4), the values lower than 1/3 are more difficult to explain. In our analysis, the exponents close 
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to 2/9=0.22 can simply be explained (in the complete condensation regime and for clean substrates) by the fact that 
the island density is not constant (as assumed in previous studies). 



V. INTERPRETATION OF EXPERIMENTAL DATA 

In principle, Figure |^ allows to determine the value of the microscopic parameters (diffusion, evaporation) if the 
saturation island density is known. The problem is : does this island density correspond to the defect concentration 
of the surface or to homogeneous nucleation? Is evaporation present in our experiments and what is the magnitude 
of r e ? 

The first question can be answered by looking at the density evolution with the flux. As already explained, if this 
leaves unaffected the island density, nucleation is occuring on defects. A similar test can be performed by changing the 
substrate temperature, but there is the nagging possibility that this changes the defect concentration on the surface. 
It is also possible to study the kinetics of island nucleation, i.e. look at the island concentration as a function of 
thickness or coverage. The presence of defects can be detected by the fact that the maximum island density is reached 
at very low coverages (typically less than 1%, see Fig. ||). One should be careful however to check that all the islands, 
even those containing a few atoms, are visible in the microscope images. 

The second question is more delicate. First, one should check whether atomic reevaporation is important. In 
principle, this can be done by measuring the sticking coefficient, i.e. the amount of matter present on the surface 
as a function of the matter brought by the beam. If possible, this measure leaves no ambiguity. Otherwise, the 
kinetics of island creation is helpful. If the saturation is reached at low thicknesses (e sa t < -5 ML), this means 
that evaporation is not important. Another way of detecting atom evaporation is by studying the evolution of the 
saturation island density with the flux : the exponent goes from 0.29 to 0.66 (Fig. ||). Suppose now that one finds 
that evaporation is indeed important : before being able to use Fig. [| one has to know the precise value of r e , 
and this is not easily achieved. For example, one could try to measure the sticking coefficient or the quantity of 
matter needed to reach saturation to obtain an estimation of the evaporation. Intuitively, the more evaporation is 
important, the more matter we need to reach the saturation density. Unfortunately, this strategy, although useful 
for growth of 2d islands Jul is not so straightforward here. The reason is that in the limit of strong evaporation 
(section [IB 2| ), one has e sa t ~ -^sat thus bringing no independent information on the parameters. The same is 



true for the sticking coefficient, which is a constant, i.e. independent of the value of r e or the normalized flux. This 
counterintuitive result can be understood by noting that in this limit, islands only grow by direct impingement of 
atoms within them. Fortunately, in many experimental situations the limit of high evaporation is not reached and we 
"benefit" from (mathematical) crossover regimes where these quantities do depend on the precise values of r e . Fig. ^ 
gives the evolutions of S sa t and e sa t for different values of r e and F in this crossover region. Then, a measurement of 
S sa t or e sa t can shed light on experimental the value of r e . For more details on interpretation of experimental data, 
we refer the reader to a review paper to appear p8j . 



VI. SUMMARY, PERSPECTIVES 



We have presented a comprehensive theoretical analysis of the growth of three dimensional structures on surfaces 
by atom evaporation. The study has been carried out by combining a simple scaling analysis, a more rigorous rate 
equations approach and computer simulations, with the main scope of helping experimentalists to analyze their data. 
The scaling analysis can give very simply the growth exponents in the "clean" substrate case, in the limiting regimes 
of high and low evaporation as well as in the crossover between these. The rate equations confirm this analysis and 
predict the growth on "dirty" substrates, i.e. surfaces containing perfect traps for adatoms. The two approaches were 
compared to Monte-Carlo computer simulations and very good agreement was found. In addition to the analytical 
predictions, computer simulations allowed to predict for the first time important growth characteristics such as 
the island size distributions and intermediate regimes which are difficult to study analytically. This is particularly 
interesting for the interpretation of experimental data from the the sticking coefficient and the saturation island 
density (see section |v|). 

We wish to thank C. Henry for helpful discussions. 
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FIG. 1. Typical island morphology generated by our model. 



FIG. 2. Universal function rescaling all our data. As predicted by Equation 2.2£ , the normalized island density N ma xX s 
follows a single curve as a function of t he ev aporation parameter (F/D) 2 ^ s Xs lu/,i (l + Xg) 2 ^ 3 . The solid curve shows the 
function predicted in the text (Equation 2.28), while the circles represent the results of the computer simulations. 



FIG. 3. Evolution of the maximum island density as a function of the flux for different evaporation times. The solid lines 
show the expected value for the exponent when evaporation is significant (2/3) while the dashed line shows the exponent in 
the complete condensation case (2/7). 



FIG. 4. Maximum island density as a function of the evaporation time for different fluxes. The number next to each symbol 
corresponds to the log(F) value for that set. The solid line shows the expected value for the exponent : 2/3 for low values of 
T e (evaporation is significant). 



FIG. 5. Exponent of the mean island radius as a function of deposition time in three deposition conditions : (a) complete 
condensation, (b) strong evaporation (c) intermediate case, (a) shows the time exponent for the radius (filled triangles) and for 
the islands (open squares). We note that the exponent of the island evolution starts at 0.33 as predicted but then slightly decreases 
to 0.26 (b) shows the averaged time exponent of the radius (open squares) as well as the real mean values of the radius for 8 runs 
(to show the fluctuations). The dashed line indicates the predicted value for the exponent (1). (c) shows the values of the radius 
exponent for two different values of the parameters, in the intermediate regime (high evaporation at the beginning of the growth, 
decreasing as islands form, see below). There seems to be a typical value for the exponent of about 0.6. In all these figures, the 
value of the local exponent is obtained by a simple derivative in the log-log plot. The precise parameters used for each graph are 
: (a) F — 10 -11 ,T e = 10 15 ,r = 1L = 1500, averaged over 5 runs; (b) F — 10 _8 ,r e = l,r = 1L = 1300, averaged over 8 runs; 
(c) squares : F = 10 _10 ,r e = 100, r = 1L = 1550, averaged over 3 runs, filled circles : F — 10 _8 ,r e = 100, r = 1L = 1350, 
averaged over 8 runs 



FIG. 6. Evolution of the island density as a function of the thickness (e = Ft) for different growth hypothesis. This figure 
shows that the same saturation density can be obtained for films grown in very different conditions. Note that the horizontal 
scale is logarithmic : therefore, nucleation on defects leads to saturation at extremely low coverages, almost impossible to observe 
experimentally. The different sets of data represent : triangles : growth with evaporation, r e = 100r and Ft = 1.2 10 -8 , circles 
: growth without evaporation (Ft = 3 10~ 10 ), and squares : growth on defects (defect concentration : 5 10 -4 per site) and 
Ft = 10~ 14 (no evaporation). 



FIG. 7. Normalized island size distributions obtained for Ft — 10 and different values of the evaporation time r e . The 
size distributions are averaged for different coverages 9 between .05 and 0.2. The solid line shows the size distribution obtained 
without evaporation. The number next to each symbol corresponds to r e /r. 



FIG. 8. Effect of the presence of defects on the island size distribution. The rescaled island size distributions are obtained 
for Ft — 10 -8 and different values of the evaporation time r e (r = 1). The size distributions were obtained for different 
coverages 6 between .05 and 0.15. Contrary to what is observed for homogeneous nucleation, the histograms do depend on the 
coverage for nucleation on defects. The solid line shows the size distribution obtained without evaporation. 



FIG. 9. Values of (a) the sticking coefficient S sa t and (b) the thickness e aa t at the saturation of island density in the total 
coalescence limit. In the limit of low island densities, S sa t is a constant. However, there are crossover regimes which depend 
on the precise r e and which are shown here. Then, from a measure of S aa t and N sat one can get an estimate for r e for the not 
too low island densities which correspond to many experimental cases. In the same spirit, (b) shows the evolution of e sa t as a 
function of N sa t in the crossover regime. The numbers correspond to the different t s /t used for the simulations and the solid 
line represents the limiting regime (see p8j for more details). 
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